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Within the atom-vacancy exchange mechanism in a nearest-neighbor interaction model we in- 
vestigate the kinetics of surface-induced ordering processes close to the (001) surface of an fee A3B- 
alloy. After a sudden quench into the ordered phase with a final temperature above the ordering 
spinodal, Tf > Tsp, the early time kinetics is dominated by a segregation front which propagates into 
the bulk with nearly constant velocity. Below the spinodal, Tf < Tsp, motion of the segregation wave 
reflects a coarsening process which appears to be slower than predicted by the Lifschitz-AUen-Cahn 
law. In addition, in the front -penetrated region lateral growth differs distinctly from perpendicular 
growth, as a result of the special structure of antiphase boundaries near the surface. Our results are 
compared with recent experiments on the subsurface ordering kinetics at CusAu (001). 
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1. INTRODUCTION 



Near a surface, kinetic processes connected with first- 
order phase transitions are generally modified in a fun- 
damental manner relative to the corresponding bulk pro- 
cess. Well-known examples are surface-induced spinodal 
decomposition and nucleation phenomena,Ei □ whose un- 
derstanding and control is of major concern in modern 
thin film and nanostructure technologies. Surface ki- 
netic effiects during ordering transitions in metallic binary 
alloysO'El acquire special interest from the scientific point 
of view as their (low-temperature) ordered phase nor- 
mally is characterized by a multicomponent order param- 
eter, whose components in general couple differently to a 
crystal surface with a particular symmetry and hence will 
display different relaxational behaviors. In such cases, 
the surface is expected to induce rich dynamical behav- 
ior of the order parameter, predominantly at short times. 

A well-studied system is the (001) surface of CusAu, 
an fcc-alloy which undergoes a first-order bulk transition 
from the disordered phase to the ordered LI2 structure 
at a temperature Tq = 663K. In the bulk, Au-atoms 
preferentially occupy one of the four simple-cubic sub- 
lattices of the underlying fcc-lattice. The ground state 
therefore is 4-fold degenerate, and there exist two types 
of antiphase boundaries separating the 4 types of en- 
ergetically equivalent equilibrium domains. Several pe- 
culiar modifications of order near the (001) surface of 
CusAu have been reported. Below the ordering temfjeL^ 
ature Tq the (OOl)-surface displays disorder-wetting.Eltil 
Moreover, above Tq, the tendency of Au-atoms to en- 
rich in the outermost (first) layer leads to an oscillatory 
segregation profile with successive Au-depletioa, aad en- 
richment in even and odd layers, respectively.tHrta This 



profile decays towards the bulk on a length scale given 
by the bulk-correlation length ^(T), which on extrapola- 
tion to lower temperatures appears to diverge at the spin- 
odal temperature Tsp ~ Tq — 30 K.E3E30 Intriguing non- 
equilibrium behavior of segregation amplitudes has been, 
detected recently by time-resolved X-ray experiments.E£l 
After a quench from an initial temperature Ti > Tq to 
a final temperature Tf < Tq the initial surface segrega- 
tion profile induces a segregation front which rapidly pro- 
gresses towards the bulk, while lateral order was found 
to develop on longer time-scales. This two-stage order- 
ing process was interpreted recently in terms of time- 
dependent Ginzburg-Landau (GL) theory.E3 Although 
restricted to one dimension and to the limit of vanishing 
thermal noise, this theory successfully yields the differ- 
ence in time-scales for perpendicular and lateral order- 
ing and identifies an anisotropy of the evolving domain 
structure within a characteristic penetration depth of the 
segregation front. 

Motivated by these studies of the subsurface ordering 
kinetics of CuaAu (001), our aim here is to gain more in- 
sight into possible surface-induced ordering scenarios in 
A3B fcc-alloys with the help of dynamical Monte Carlo 
simulations. As a minimal model with respect to equilib- 
rium properties wc choose a lattice with nearest-neighbor 
interactions among A and _B-atoms, known to account 
reasonably for the bulk phase diagram of the Cu — Au 
system. Our dynamic algorithm relies on the vacancy 
mechanism, where elementary atomic moves consist in 
anpfipoihange of (^)Cu- or (i3)Au-atoms with vacancies 
T/.ETlD This kind of ABV-mode\, which has been used 
before by Frontera et al. in arLiiivestigation of the bulk 
ordering kinetics in A^B alloys Ba is physically more real- 
istic than the conventional direct Ai3-exchange kinetics 
and will simultaneously allow us to scale the simulated 
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surface-induced growth rates with the self-difFusion co- 
efficient of A- and _B-atoms. 

The main results of our simulations are: (i) the oc- 
currence of two distinctly different modes of progression 
of the segregation front, depending on the final temper- 
ature Tf to lie significantly above the spinodal temper- 
ature, Tf > Tsp (quench into the metastable regime) or 
below, Tf < Tap (quench into the unstable regime), (ii) 
an early-time linear (constant velocity) growth behavior 
of the segregation front in the case Tgp <Tf <To and a 
characteristic time t*, where this linear growth deceler- 
ates. Below Tgp it changes over to a growth roughly pro- 
portional to i^/^ up to our longest computing times, (iii) 
A different behavior of lateral growth in the near-surface 
region. In connection with this we give a detailed de- 
scription of the persistent anisotropy in the near-surface 
domain-structure. 

The paper is organized as follows. Section 2 defines 
our model; its bulk properties and equilibrium surface 
properties are tested in Section 3. The subsequent Sec- 
tion 4 discusses in detail the processes of perpendicular 
and lateral ordering. Our findings have strong relevance 
to the temperature quench experiments on CusAu de- 
scribed above, although the present model based on spa- 
tially uniform parameters for nearest-neighbor interac- 
tions cannot be expected to work quantitatively. In fact 
it implies an enhancement of static lateral order near the 
surface as compared to the bulk, an effect not borne out 
in the real material CusAu. Enhanced lateral order at 
the (OOl)-surfaceJiowever, seems to be observable in the 
fcc-alloy CusPd.Ea Concluding remarks are presented in 
Section 5. 



2. MODEL AND SIMULATION METHOD 

Consider a lattice oi N x N x M fcc-cells with periodic 
boundary conditions in the x- and y-direction. Along 
the 2;-axis our system is bound between two free surfaces. 
Atomic layers normal to the z-direction with distance a 
are labeled by n with 1 < n < 2M + 1. Each lattice site 
i can be occupied either by an A-atom, or a _B-atom, or 
a vacancy so that the corresponding occupation numbers 
cf, cf and satisfy cf + cf + cY = 1. In our simphfied 
alloy model only pairwise nearest-neighbor interactions 
are taken into account. The configurational Hamiltonian 
then reads 

H^Y. [y^Actcf -f VBBcfcf + Vab [ctcf + cfcf)] , 

(i,J> (1) 

where the summation is over nearest-neighbor pairs.0 In 
view ofjthe very small concentration of vacancies in real 
alloystll it is natural to assume c( <^ 1 such that macro- 
scopic equilibrium properties of the alloy phases will re- 
main essentially unaffected by the vacancies. For equi- 
librium considerations one can then ignore the vacancies 



and introduce spin variables Si = 2cf — 1 = ±1 by which 
equation (Q) is mapped onto the spin-1/2 Ising model, 

Hi = -jJ2 -hJ2s^-hlJ2 (2) 

with 

J=-\{Vaa + Vbb-2Vab) (3) 

h = SiVBB - Vaa) (4) 

The last summation in is restricted to spins in the sur- 
face layers n = 1 and n = 2M + 1, which are subjected 
to the surface field 

hi = Vaa ~ Vbb (5) 

Bulk equilibrium properties at fixed composition only 
depend on the parameter J. Regarding the Cu(yl)- 
Au{B) system, it is natural to take Vbb > because 
of the larger size of Au-atoms relative to the Cu-atoms, 
whereas Vaa < 0, Vab < such that bulk interactions 
will favor "anti-ferromagnetic" ordering J < 0. Near the 
stoichiometric composition A^B the emerging bulk LI2- 
structure is described by a 4-component order parame- 
ter, one conserved concentration variable ipo and three 
non-conserved components i/'ij "02, and "03 • These are de- 
fined via the differences ttIq = {cf) — (cf ); i G a; of 
averaged A— and _B— occupations of the four equivalent 
simple-cubit. sublattices a = 1, . . .4 building the fee lat- 
tice, namclytl 

■00 — TTll + 7712 + ITilj, + 7714 
01 = TTll - 777,2 - 7773 + 7774 
W2 = 777.1 — 7772 + 7773 ^ '^74 
03 = 777.1 + 7772 — 7773 ~ '^74 

Each non-conserved component describes a layerwise 
modulation of the _B-atom concentration along one of 
the three cubic axes. The ground state shows 4- 
fold degeneracy, where the four types of equivalent 
ordered domains have components (0i, 02, 03)/'0 — 
(-1,1,1); (1,-1,1); (1,1,-1) and (-1^-1,-1). In the 
ground state (T = 0) the amplitude takes the value 
^ = 2. 

The above conditions for the interaction parameters 
imply a negative surface field hi < 0, favoring an en- 
richment of B-atoms in the outermost layers 77 = 1 and 
n — 2M + 1. This is to be expected in view ofJJie BB- 
repulsion. As shown by Schweika and Landauu on the 
basis of the Hamiltonian and single-spin flip dynam- 
ics, the value hi/\J\ = —4 is consistent with the experi- 
mentally measured concentration of Au-atoms of about 
50% in the surfajce layer of Cu3Au(001) near the ordering 
temperature ToS3 

In our present work, all investigations will be based on 
the atom-vacancy exchange mechanism, where the en- 
ergetics are given by the complete hamiltonian (1^). For 
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the final choice of interaction parameters see the next 
section. To affect an elementary atomic move in the 
Monte Carlo runs we randomly select a nearest-neighbor 
vacancy-atom pair and calculate the vacancy-atom ex- 
change probability by using the Metropolis algorithm. 
Throughout most of our calculations the mean vacancy 
density is taken as ~ 6.1 • 10~^, which amounts to 128 
vacancies in a system of typical size TV = 64; M = 128 
with m'^M sites. One Monte Carlo step (MCS) consists 
of AN'^M attempted moves and is independent of c . 
Clearly, in an actual material the vacancy concentration 
generally will depend on temperature and the amount of 
A- and i3-atoms. Our kinetic model with fixed there- 
fore does not allow us to compare time-scales of vacancy- 
mediated processes taking place at different temperatures 
and composition. However, our subsequent studies of 
the ordering kinetics focuses on the exact stoichiometric 
composition A^B (c^ = 0.25) and a narrow temperature 
range near Tq where changes of are expected to be of 
minor importance. Local order parameters to be used in 
the analysis of our simulations are defined in terms of av- 
erages of over one fcc-cell. Order parameter profiles 
and structure factors are obtained by averaging over 10 
independent simulation runs. 



3. BULK PROPERTIES AND SURFACE 
SEGREGATION 

In this section we discuss some elementary properties 
of our model, part of which are already known in the lit- 
erature but they are needed here as a test of our kinetic 
algorithm and as a basis for the subsequent investigations 
of surface-induced kinetic properties. Fig. |l|a shows the 
phase diagram in the vicinity of ~ 0.25 as deduced 
from simulations of a system with size M ~ N ~ 64 
and periodic boundary conditions in all directions. To 
determine the ordering temperature at stoichiometry 
[c^ = 0.25), we analyzed the decay of a perfectly or- 
dered initial state upon thermalization at different tem- 
peratures. This yields fcs To ~ 1.83| J| . in goo d agreement 
with literature data (fcsTo ~ 1.85| J|).EaL3S From a sim- 
ilar procedure, supplemented by computations of concen- 
tration (?/'o)-histograms, we deduced estimates for the 
boundaries of the two-phase regions, which reasonably 
reflect the bultnphase diagram of the Cu — Au system 
near = 0.25.eJ 

Next we proceed to equilibrium properties connected 
with the (OOl)-surface. From our algorithm defined in 
Section 2 we computed equilibrated _B-atom concentra- 
tions in the surface layers as a function of /ii, keeping J 
fixed. As before, hi and J are defined by (J^) and (|^). 
Good agreeroent was found with the work by Schweika 
and Landau.o We consider this both as a positive test 
of our algorithm and as an indication that the small 
amount of vacancies used has practically no influence on 
these static results. In what follows we shall therefore 
choose the interaction parameters of our CusAu-model 



such that hi = — 4|J|.e3 While J and hi are now de- 
termined via the bulk ordering temperature Tg and the 
Au-segregation amplitude in the surface layer one more 
interaction parameter needs to be specified. Our final 
choice is Vbb — —Vaa — —Vab > with |J| = Vbb/'^-- 
This can be shown-ip-jmply a weak effective attraction 
between vacancies. tJ'EZl 

Furthermore we remark that the present model dis- 
plays a surface transition at a critical temperature Tcs > 
To, which we estimate as T^s ^ 1-12 To for = 0.25. Be- 
low Tcs the outermost layers with nearly equal amounts of 
A- and i3-atoms develop a lateral (l,l)-superstructure. 
With respect to CuaAu, the existence of a surface critical 
temperature clearly is an artefact of the present model, 
connected with the fact that fee lattice models based 
on spatially uniform short range interacticps are unable 
to account for surface-induced disordering .b3 The (001)- 
surface of the Ll2-structure of CusPd, on the other hand, 
appears to favor lateral ordering.c3 
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FIG. 1. (a) Phase diagram of the ^_B-binary alloy model 
in the vicinity of = 0.25. Data points from simulations 
are connected by smooth lines. Full line (■) and dashed 
line (•): boundaries of the disordered phase and the ordered 
I/l 2 -structure, with the two-phase coexistence region in be- 
tween. Dotted line (a): ordering spinodal. 
(b) Plot of versus temperature at — 0.25, illustrating 
the determination of the spinodal temperature. 
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Moreover we estimate the spijiodal temperature in a 
manner as done experimentallyJ13 In the disordered phase 
B-atoms are known to show an osciUatory segregation 
profile near the (001) surface, with enrichment and de- 
pletion in odd and even layers, respectively. When atoms 
in the outermost layers are identified with sublattices 
a = 3 and 4, the segregation amplitude is given by 
the order parameter component ■03. At the surface, its 
magnitude is determined by the surface field hi. Below 
the surface it decays towards the bulk on aJength scale 
given by the bulk correlation length ^(T).E3 Extrapola- 
tion of data for ^ (T) to lower temperatures according to 
^(T) ~ (T — Tsp)"-*^/^ allows one to estimate the spin- 
odal temperature Tsp. Experimentally, Tq ~ 663 K and 
Tsp ~ 30 K near the stoichiometric composition of 
Cu3 Au.ll3 Application of this procedure to our simula- 
tion data is illustrated in Fig. and yields the ordering 
spinodal displayed in Fig. |]a. At = 0.25 we find 
Tsp/To ~ 0.967 ± 0.003.EI This is somewhat larger than 
the experimental value but reasonably close in view of 
the simplicity of our model. Simultaneously we get (see 
Fig. |]b) S.i'fe^ — 6a, which again agrees fairly well with 
experiment J13 
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FIG. 2. Bulk domain patterns (c — 0.25) at a time 
t = 7 ■ 10^ MCS after a quench from random initial condi- 
tions. Grey scales reflect values of the local order parameters 
between extremal values 4 and —4. 

(a) i/^i-pattern at a final temperature T/ — 0.972 To > Tsp. 

(b) -(d) patterns of V"!; V'2, V's at Tf = 0.963, making evident 
the two types of domain walls and the anisotropy of the evolv- 
ing structures. 



For metallic alloys where short range interactions pre- 
vail, it is well known that tie spinodal and the nucleation 
regime connect smoothly.cl The qualitative significance 
of the spinodal temperature, however, becomes apparent 
from Fig. |^a and b. There we compare bulk patterns of 
the order parameter component ipi at a time t = 7 ■ 10'^ 
MCS after a sudden quench from a purely random initial 
state to final temperatures T/ ~ 0.972 Tq > T^p (Fig. |a) 
and Tf = 0.963To< Tsp (Fig. |b), respectively Snap- 
shots were taken in a section parallel to the {xy)-p\a,ne. 
In the first case only small ordered domains have nucle- 
ated in an essentially disordered background (Fig. |^). 
By contrast in Fig. ^b we observe a continuous domain 
pattern typical of spinodal ordering. At that temperature 
{Tf — 0.963 To) the same type of patterns is seen also at 
shorter observation times. For a complete characteriza- 
tion of the evolving domains in Fig. ^b we need in addi- 
tion the ^p2~ and ■(/'3-patterns displayed in Figs, and 
From the grey-scales in Figs, ^s-^jd one can easily see 
that inside a domain we have ipiip2''p3 < 0- Furthemore, 
the well-known two types of antiphase boundariesLj are 
recovered. For example, boundaries where V'l and ip^ si- 
multaneously change sign in the y-direction are the "low- 
energy" (type-1) domain walls and in fact cost no energy 
in our nearest-neighbor model. These walls are almost 
completely flat and turn out to be rather stable. On the 
other hand, boundaries where ipi and ip3 simultaneously 
change sign in the a:-direction are "high-energy" (type- 
2) walls where curvature-driven coarsening is much more 
effective than for type-1 walls. Bulk coarsening effects at 
quench temperatures lower than those considei^d here 
have been simulated in detail by Frontera et al.EEl These 
authors obtained coarsening exponents n smaller than 
that in the conventional Lifschitz-AUen-Cahn coarsen- 
ing law, where n = 1/2, and interpreted their findings in 
terms of the presence of those stable low-energy walls. 
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FIG. 3. Temperature-dependent diffusion constants of A- 
and B-atoms in an Arrhenius representation. 

In addition we extract tracer diffusion coefficients from 
our model, restricting ourselves again to = 0.25 and 
to temperatures above but close to Tq. Diffusion coef- 
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ficients Da and Db for A- and i?-atoms are obtained 
in the usual way from time-dependent mean-square dis- 
placements and are plotted in Fig. ^ in an Arrhenius rep- 
resentation. Since in our model Vaa = Vab, the energy 
remains essentially unchanged by a jump of an A-atom so 
that Da is found to be nearly temperature-independent. 
On the other hand, because of the B-i3-repulsion the B- 
atoms will migrate more slowly than the ^-atoms, and 
Db decreases upon cooling. As indicated already, in a 
real system will grow with temperature. Therefore 
the individual experimental curves log Da or log Db ver- 
sus 1/r will be steeper than those in Fig. ||. However, 
the ratio Da/Db roughly agrees with experimerital data 
for CusAu, where Da/Db ^ 1.45 at T/Tq ~ 1.7^ while 
we find Da/Db — 2 at the same T/To-ratio. 

Let us close this section by a brief remark on the dis- 
tribution of vacancies in our system with free surfaces. 
First of all, vacancies tend to enrich in the outermost 
layers. For example at T/Tq = 0.972 we find a surface 
vacancy concentration of about c^^^f^^^ ~ 2.5-10^^, much 
larger than the overall concentration. The fact that va- 
cancies are expelled from the bulk suggests that in fact 
we are dealing with a model for a stable solid which, 
when regarding the complete ternary {ABV)-system, 
would be phase-separated from the "vacuum" (vacancy- 
rich) phase. Furthermore, when analyzing the above- 
mentioned layered structure along the z-axis caused by 
surface segregation, we find that vacancies deplete in the 
ABAayers and enrich in the v4-layers. The reason is 
that the BB repulsive energy in ABAayers becomes min- 
imized when vacancies avoid i?-sites. The resulting os- 
cillations in the vacancy-density along the z-axis have 
an amplitude of about 30 percent relative to the average 
vacancy density. 



4. SURFACE-INDUCED ORDERING KINETICS 

Now we turn to the question of how order evolves in 
a system with free surfaces. Quench conditions are cho- 
sen as in Section 3 with final temperatures below Tq but 
in the vicinity of T^p {0.9<Tf /Tq^I). First, within few 
MC-steps i3-atoms enrich in the surface layer to a con- 
centration of about 50% and deplete in the adjacent layer, 
thereby establishing local equilibrium as enforced by the 
surface field hi. This implies Va — 1-85 to about 1.9 at 
2; = in the temperature interval considered. This pro- 
cess in turn induces moves of _B-atoms from even to odd 
layers within an increasing depth from the surface, thus 
forming a '03-front which penetrates into the bulk. As an 
example, we show in Fig. ^ the ■i/'a-distribution in a sec- 
tion parallel to the {y, z)-plane at Tf = 0.972 Tq > Tgp 
for three different times after the quench. The patterns 
clearly demonstrate a systematic progression of an or- 
dered state with fixed (positive) sign of i/'s- Taking aver- 
ages over sheets parallel to the surface, we obtain profiles 
ip3{z, t) which are plotted in Fig. |[ 
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FIG. 4. Vs-patterns at Tf = 0.972 Tq for three different 
times after the quench, illustrating propagation of a sur- 
face-induced ordering wave. 

At t = 10"* a plateau region for smaller z has evolved 
with V'3~values in the vicinity of the equilibrium order 
parameter ip ~ 1.79 at that temperature, before the 
profile decays to zero for larger z. Evidently, the thick- 
ness of the interface between the partially ordered state 
and the bulk (not penetrated by the tps-iront) increases 
with time, an effect expected to arise both from growing 
intrinsic capillary wave fluctuationtca and from bulk fluc- 
tuations in all three order parameter components with 
the character as in Fig. ^ja, to be "incorporated" into the 
front during propagation. 
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FIG. 5. Average profiles tp3{z,t) for several times after the 
quench. 

The profiles ip^lzjt) allow us to define a time- 
dependent penetration depth Z3{t) by setting 
^J3{z3{t),t) = ^ip3{0,t). Plots of Z3{t) for a series of 
final temperatures Tf are shown in Fig. ^ for a system 
large enough in the z-direction {N — 64, M = 128) to 
avoid the ^/^a-fronts induced by both surfaces to over- 
lap. (In fact, data in Fig. ^ and |^ contain an average 
over those two fronts.) As expected, after a quench not 
passing the phase boundary 2:3 (i) rapidly saturates. For 
example, at Tf = 1.006 Tq the final penetration depth is 
about 6 atomic layers, which is of the order of the bulk 
correlation length ^ in the disordered phase at that tem- 
perature. On the other hand, when Tsp < T/ < Tq, we 
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see that after a short transient Z3(t) grows more or less 
Hnearly with time within an extended time interval. This 
allows us to introduce a penetration velocity v{T) which 
increases upon lowering the temperature. For example, 
at Tf — 0.972 To we find that 23 (i) grows linearly up 
to at least 10'* MCS, where it reaches about 70 atomic 
layers, see also Fig. ^. In determining v(T), shown in 
Fig. 0, care has to be taken with respect to the lateral 
size N of the system. For small N the system is likely to 
develop single (closed) domains in the lateral direction 
within the front region. Increased order of this type will 
speed up the penetration process. This iV-dependence of 
V is exemplified in Fig. |^ for Tf — 0.972 Tq. Data up to 
t — 10* MCS, however, seem to converge sufficiently at 
our largest N — 128. Moreover, using again N — 64, we 
have also tested the influence of the vacancy concentra- 
tion on V. At Tf = 0.972 Tq, results with = 3.05 ■ 10"^ 
(64 vacancies) are indistinguishable from those shown in 
the figure. Hence the vacancy system appears to be in 
the dilute limit also with respect to growth processes of 
this type. Note that for CusAu the-experimental vacancy 
concentration at Tq is even lower.E2l 

As seen from Fig. ^ when temperature is decreased 
from To, the velocity v{t) first increases with upward cur- 
vature, but decelerates when T « Tgp. Sufficiently below 
Tsp, linear growth can no longer be identified and w(T) 
becomes undefined, see the case Tf — 0.911 Tq in Fig. m. 
Some estimates based on mean-field arguments concern- 
ing the order of magnitude of w(T) are presented in the 
Appendix. 

When temperature falls near or below Tf ~ 0.963 Tq, 
a new regime of slower growth beyond some crossover 
time t* becomes detectable within the time window of 
Fig. ^. t* decreases with temperature, so that the regime 
of linear growth shrinks. Notice in addition the non- 
monotonous behavior of 2:3 (t) as a function of tempera- 
ture for fixed observation time t. As indicated already, 
the strong shrinkage of the linear regime, where t* tends 
to merge with the short time transient, occurs at temper- 
atures slightly below but very near the spinodal temper- 
ature Tsp ~ 0.967 discussed in Section 3 on the basis of 
static considerations. Regarding these surface phenom- 
ena we conclude that the spinodal temperature as esti- 
mated via the static segregation profiles separates two 
distinct dynamical regimes, roughly characterized by the 
presence or the absence of surface-induced linear (con- 
stant velocity) growth in an intermediate time domain. 

The slowing down of the front motion beyond the lin- 
ear regime, i.e. for t > t* , clearly originates from its 
competition with nearly saturated domains that have al- 
ready nucleated in the bulk. Crudely speaking, t* will be 
determined by the interplay of two dynamic quantities, 
the penetration velocity v{T) and the nucleation rate of 
ordered domains. Moreover, unstable growth of bulk do- 
mains at temperatures Tf < Tsp totally suppresses the 
linear regime. To investigate this situation more closely, 
we carried out even longer runs up to t = 2 ■ lO'' for the 
case Tf — 0.911 Tq. Results at that temperature for z^{t) 



in a log-log representation are shown in Fig. The data 
in the range t^2 ■ IQ^ MCS appear to follow a power-law 
z^{t) ~ t" with n ~ 0.27 ± 0.02. 
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FIG. 6. (a) Time-dependent penetration depth 23 (t) of the 
i/ia-front (see Fig. ^) for various final temperatures. 

(b) Dependence of 23 (t) on the lateral system size A''. 

(c) Time-dependent penetration depth in a double-logarith- 
mic representation for three different temperatures. The 
growth behavior for large times (up to t = 2 • 10 MCS for the 
two lower temperatures) is compared with a t^^^'-law (straight 
line). At Tf = 0.963 Tq a comparatively short time domain 
of linear growth appears (see the steeper straight line) up to 
the crossover time t* . 
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The appearance of such small growth exponent is con- 
firmed by additional simulations at a temperature Tf — 
0.759 To considerably below Tgp, also shown in Fig. 
These data indicate in addition that the onset of this 
type of power law is moved to earlier times when tem- 
perature is lowered. Conversely, data for Tf — 0.963 Tq 
display a crossover to a slow temporal growth compatible 
with n ~ 1/4 at a larger crossover time, indicated as t* in 
Fig. ||c. Since the penetration depth up to i = 2-10^ stays 
significantly smaller than the system size (in this case 
N = M = 128), we expect that finite size effects do not 
play an essential role in estimating that exponent n. Pre- 
liminary results for the structure factor of a bulk system 
seem to confirm that value of n, although no definitive 
conclusions concerning the asymptotic long-time behav- 
ior can be drawn from these studies so far. Nevertheless 
it is interesting to compare our findi ngs i n Fig. ^ with 
recent work by Castan and Lindgard.r'irH These authors 
studied coarsening processes in a 2-dimensional system 
involving both curved and flat (zero curvature) domain 
walls. The latter were essentially immobile and could 
be removed only by the progression of an intervening 
curved wall. Such a situation basically is encountered 
also in the present model, where type-1 antiphase bound- 
aries turn out to be essentially immobile. This stability 
of fiat walls induces an overall growth exponent which 
is smaller thajXuthe classical Lifschitz-AUen-Cahn expo- 
nent n — 1/2. oe3 Castan in fact gave arguments in favor 
of an exponent n = 1/4. A similar conckision with re- 
spect to n was reached by Deymier et al.o from studies 
of a different model with two types of contacts between 
domains. 



Our main objective, however, is the near-surface region 
z < Z3{t) and consequences of the s urfa ce-induced segre- 
gation wave. As emphasized before,t§Ej for z < z^lt) the 
sign of ip3 is essentially fixed to ^/'a > so that only two 
types of domains can appear: {'ip 1,1^2,1^3)/ 4' = (1, —1, 1) 
and (—1,1,1). These correspond to the two ways B- 
atoms can form a (1, l)-superstructure on the square 
lattice building an odd (AB-type) layer. Walls between 
those domains encountered when going in the lateral di- 
rection are exclusively of type 2, whereas type-1 walls are 
parallel to the surface. To analyze this situation further, 
we have enlarged the system in the lateral direction to 
N — 256, but restricted ourselves to M = 8 so that tps 
is nearly uniform along the z-axis. Patterns of all three 
order parameter components were extracted from atomic 
layers 7 and 8 midway between the two surfaces. An ex- 
ample with Tf = 0.911 To and t = 10* MCS is shown in 
Fig. H In fact, domain walls in Figs. ||a and ^3 with si- 
multaneous sign changes of tpi and ip2 are always curved, 
indicative of type-2 walls, and the domain structure ap- 
pears completely isotropic, in contrast to the anisotropic 
shape of bulk domains (cf. Fig. ||). Fig. ^ depicts 
which is nearly constant. 
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FIG. 7. Penetration velocity v{T) versus temperature. 

Next we study lateral ordering within layers parallel 
to the {xy)-p\a,ne at a depth z from the surface. This 
process proceeds in different ways depending on whether 
z < 2:3 (t) or z > 2:3 (t). When analyzing the layer z = M 
midway between the two surfaces, representative for the 
latter case, we recover the characteristics of bulk order- 
ing as described briefiy in Section 3, see also Fig. ||. 



FIG. 8. Lateral domain patterns for a) and b) i/)2 in the 
near-surface penetrated region with nearly uniform ip^ (see 
c)) at Tf = 0.911 To, t = 10* MCS. Notice the absence of 
type-1 walls (see text). 

The time-dependent average size of domains in the 
(xy)-section of Fig. || is analyzed in Fig. ^, where we 
present results for the first moment 



fc(j(z, t) 



Efcii fc||S'a(fc||,2,i) 



(7) 



Efc|| Sa{kii,z,t) 

of the time-dependent lateral structure factor, defined by 

^„(/c||,z,f) = (|^a(fc,|,z,t)p) (8) 

where a = 1 to 3 and ipaikw, z,t) is the lateral Fourier 
transform of the local order parameter il}a{r,t). Note 
that vacancy depletion in j4i3-type layers and enrich- 
ment in A-type layers (see Section 3) expands the over- 
all time-scale for lateral ordering. While the behavior 
of k'^{t) reflects a rapid transition to a near-uniform ?/'3- 
distribution (see Fig. ^), the data points for ka{t) with 
a = 1,2 in Fig. ^ seem to approach a growth behavior 
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k~^{t) described by the Lifschitz-AUen-Cahn ex- 

ponent n = 1/2. In fact, our results indicate curvature 
driven coarsening in the lateral directions, distinctly dif- 
ferent from the perpendicular growth displayed in Fig. 




FIG. 9. First moments ka{t) of lateral structure factors 
5a(fc||,z,t) in the same (a;, 2y)-section as in Fig. ^, for 
Tf = 0.911 To. The behavior of ki{t) = k2{t) at long times is 
compared with the Lifschitz-AUen-Cahn prediction (straight 
line) . 



Moreover, our present studies exemplify different or- 
dering scenarios that may emerge more generally in sys- 
tems whose bulk structure requires a multicomponent or- 
der parameter, and where the surface favors some kind 
of partial ordering. Concerning domain growth in our 
model, we have to distinguish between perpendicular 
and lateral growth modes, between a near-surface par- 
tially ordered region (z < 23 (i)) and the bulk region 
(z > z^it)), and between quench temperatures Tf above 
and below the spinodal temperature T^p. At least three 
distinct coarsening schemes at intermediate or long times 
were observed. If Tsp < 1/ < To, perpendicular order ini- 
tially grows linearly with t such that Z3 (t) ~ vt for t < t* , 
while for t > t* our data indicate a temporal regime 
where 23(i) ~ t^^^. On the other hand, for Tf < Tsp 
the regime of linear perpendicular growth is suppressed. 
Analysis of the typical domain structure in the lateral 
direction shows that in the region z < zsit) it is isotropic 
and coarsens according to k^^{t) ^ t^^^ at long times 
(a = 1; 2) while for z > Z3{t) one recovers bulk behavior. 

Besides some open questions addressed already above 
we like to point out that non-stoichiometric alloys should 
in principle display even more rich and interesting be- 
haviors, in particular in cases where ordering and phase 
separation occur simultaneously.CJ 



5. SUMMARY AND CONCLUDING REMARKS 
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Surface-induced ordering processes in an fee alloy 
model were studied on the basis of the atom-vacancy 
exchange mechanism. One essential result of our sim- 
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induced at the surface and propagating into the bulk im- 
mediately after a temperature quench across the ordering 
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CU3AU (001) observed in recent X-ray experiments. In 
the region z < z^it) covered by the segregation wave, the 
present model with nearest-neighbor interactions only 
shows an increased tendency to build up lateral order. 
The experimental verification of this effect may depend 
on the alloy material under consideration. 

A second important conclusion confirmed by our simu- 
lations is the existence of an anisotropic domain structure 
in the near-surface region z < 2:3 (i). Only two types of 
domains and, for z fixed, only high-energy domain walls 
appear, as opposed to the well-known domain structure 
in the bulk with 4 types of domains and both high- and 
low-energy walls in either direction. Since in CU3AU z^{t) 
can become larger than 200 M3 or about 50 atomic lay- 
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ers and overlapping '(/'s-fronts induced by both surfaces. 
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APPENDIX 

We attempt here to give some qualitative estimates 
concerning the magnitude of the penetration velocity 
v(T), see Section 4, from simplified mean-field considera- 
tions. Our starting point are time-dependent Ginzburg- 
Landau equations, based oft, the free energy density 
/(V"!; "02, "03) as given by LaiEZi for the CU3AU structure, 



E 




(9) 



with r = rQ{T — Tsp);rQ > 0;u > 0;\v\ < u;w > 0. 
The ordering temperature Tq is determined by r{To) = 
(2/9)w^/ (3w + v) and the order parameter at Tq has the 
magnitude j-i/'al — "0 with ip = (2/3)w(3m -|- v). Regard- 
ing the penetration of the segregation wave as a one- 
dimensional process along the z-axis (perpendicular to 
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the surface), we assume an equation of motion for V-'3(^i t) 
of the form 

a.^3^-r(^-x0^') (10) 

with a kinetic coefficient F and K > 0. Two limiting 
cases concerning the couphng of '03 with tpi and ip2 can 
be considered such that (^ becomes a closed equation 
for '03 only. Thereby it turns out that the simulated 
velocity falls between these two limits. First, if the pa- 
rameters are such that the evolution of lateral order, ex- 
pressed by Tpi and '02, considerably lags behind the pro- 
gression of the '03-front, one may set in ( p^ 0i = '02 = 0. 
In this case, the form (^) of the free energy density al- 
lows a nontrivial solution for a propagating front only 
at temperatures T < T^p, with a velocity of propagation 
V oc {Tsp — r)^/^.lll Correspondingly, we obtain v = 
for T > Tsp, contrary to Fig. ^. However, simulations 
of our nearest-neighbor interaction model indicate that 
lateral order at least on small length scales evolves al- 
most simultaneously with the front motion in the per- 
pendicular direction. This suggests to consider a sec- 
ond, in fact opposite limit, namely that '0i and '02 in- 
stantaneously relax towards the "local equilibrium con- 
dition" = —'02 = ±'03. This yields the approxima- 
tion wipiip2 — —wipi, to be used in (p^). Introducing 
'y' = 03/'0 and the "potential" 

Viip) = -//(3r(To)V?) = 

= i(l-r/r(ro))v.^- 1^2(^-1)2 (11) 

we obtain 

= Tr{To){dV/d^ + (e(To)) V") (12) 

where — {K/r(To)y/^ denotes the correlation length 
at Tq. The ansatz ip(z,t) = rj{z — vt) leads to an or- 
dinary differential equation for s^z). The velocity v is 
determined in the usual mannero by requiring that ap- 
propriate boundary conditions can be fulfilled. These are 
r](z) -> ?7niax; ri'{z) ^ as z ^ -cx), where 77„iax > cor- 
responds to the absolute maximum of V{7f) and ri{z) 0; 
77' (z) ^ as z — > cxD. Numerical solution of this prob- 
lem shows that, at Tq, v{T) starts to grow linearly in 
To - T and that v{Tsp) ~ 1.4Fr(ro)Co- To estimate F, 
one can formulate mean-field arguments for the diffusion 
of the slower atomic species, in this case the i?-atoms, 
along the z-axis, which yields Db/cl^ — 2Fr(To). Us- 
ing ^0 — 6a, this gives the relation v{Tsp) ~ A{DB/a), 
independent of the Monte Carlo time unit. From the 
simulated value Db ^ 0.04a2/MCS at Tq, see Fig. |, we 
obtain v{Tsp) ~ 0.16a/MCS. Comparison with Fig. g 
shows that these simple arguments overestimate the ve- 
locity at Tgp by about one order of magnitude. This is 
not surprising in view of the simplifications made, espe- 
cially in view of the neglect of lateral fluctuations. These 
may strongly reduce the penetration velocity, an effect 



already inferred from the iV-dependence of the velocity 
displayed in Fig. ||b. 
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